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ABSTRACT 

We perform numerical calculations of the expected transit timing variations (TTVs) induced on a 
Hot- Jupiter by an Earth-mass perturber. Motivated by the recent discoveries of retrograde transiting 
planets, we concentrate on an investigation of the effect of varying relative planetary inclinations, 
up to and including completely retrograde systems. We find that planets in low order (E.g. 2:1) 
mean- motion resonances (MMRs) retain approximately constant TTV amplitudes for < °i < 170°, 
only reducing in amplitude for i > 170°. Systems in higher order MMRs (E.g. 5:1) increase in TTV 
amplitude as inclinations increase towards 45°, becoming approximately constant for 45 < i < 135, 
and then declining for i > 135°. Planets away from resonance slowly decrease in TTV amplitude 
as inclinations increase from to 180, where-as planets adjacent to resonances can exhibit a huge 
range of variability in TTV amplitude as a function of both eccentricity and inclination. For highly 
retrograde systems (135° < i < 180°), TTV signals will be undetectable across almost the entirety 
of parameter space, with the exceptions occurring when the perturber has high eccentricity or is very 
p j ■ close to a MMR. This high inclination decrease in TTV amplitude (on and away from resonance) is 

important for the analysis of the known retrograde and multi-planet transiting systems, as inclination 
• \ effects need to be considered if TTVs are to be used to exclude the presence of any putative planetary 

4^ ■ companions: absence of evidence is not evidence of absence. 

Subject headings: celestial mechanics — planetary systems — methods: numerical 
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1. INTRODUCTION 

Of the known extra- Solar planets, more than 60 tran sit the host star. Of these sys tems, at least fo ur show evid ence 
for an ext ernal companion (G J 436. lManess et al.ll2007l : HAT-P-13. l"Bakos et al-ll2Q09b HAT-P-7. rWinn et al.ll2009l and 
CoRoT-7 lQueloz ^tIHl200^ . 

If the transiting planet were the only planet in the system, then the period between each successive transit would 
be constant (neglecting complicating effects such as general relativity, stellar oblateness and tides). The presence of 
additional planets in the system (which themselves may or may not transit the star) can cause perturbations to the 
orbit of the transiting planet, leading to detectable transit timing variations (TTVs) of the known transiting planet 



orbit ot trie transiting planet, leading to detectable transit timing variations (1 1 VsJ ot tne known transiting planet 
(|Miralda-Escudell2002l : iHolman and Murravll2005[ lAgol et al.ll2005li. Oth er studies have extended t he TTV method to 
demo nstrate the feasibility of detecting planetary moons (|Kippingll2009d I&l) and Trojan companions (jFord and Holmanl 
\ l2007f ) in extra-solar planetary systems. 

However, practical use of TTVs as a detection tool requires the solution of the difficult "inverse problem" , i.e. given 
a particular TTV profile, can one reconstruct (or at least restrict) the mass and orbit of the unseen perturber? This 
prob lem is non- trivial, as numerous different perturber mass-orbit configurations can lead to degenerate TTV solutions 
(e.g. lNesvornv and Morbidellill20M ). 
^ ■ Several investigations have considered (and ruled out) planets more massiv e than the Earth in certa i n orbits / regions 
of parameter space close to many of the known transiti ng planets (e.g. ISteffen and Agoll [20051; lAgol and Steffenl 
, & ! 120071 : lAlonso et all 120081 : iBean and Seifahrtll2008l ). while iNesvornv and Morbidellil (j2Q08h and INesvornvl f2009l ) have 
developed an approximate analytic method to try and tackle the inverse problem in a more general manner. 

However, most previous investigations have focused on analyzing the effect of prograde, c oplanar companions. Even 
where explicit investigations of inclination effects have been conducted (e.g. INesvornvl [2009). the investigations have 
been restricted to prograde cases, examining only relative inclinations significantly lower than 90°. 

There are several reasons to consider TTVs of inclined systems. Observationally, it is possible to measure the sky- 
projected angle between the spin vector of a star, and the orbital angular momentum vector of a planet transiting that 
star via the Rossiter-McLaughlin effect. Measurements over the past year h ave revealed t hat a num ber of the known 
transiting planets are strongly inclined, or e ven retrograde (e.g. HAT-P-7b. [Narita et al.ll2009l and lWinn et al.l l2009. 
as well as WASP- 17b. I Anderson eralll2009f ). 

Although the mechanism(s) driving the creation of such retro grade planets remains unclear, models used to ex- 
plain high eccentricity exo-planets via planet-planet scattering (jRasio and Ford 1996; Weidenschilling and Marzari 
1996; Lin and Ida 1997; Levison et al. 1998; Papaloizou et al. 2001; Moorhead and Adams 2005; Chatteriee et al. 2008; 
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Juric and Tr emaine 2008: iFord and Rasioll2008[ ) and/or Kozai oscillations (E.g. iKozail 119621 : iTakeda and Rasiol 12005: 
Nagasawa et al.l l2Q08h can naturally excite large orbital inclinations (although the quantitative details of the distri- 
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butio ns can differ greatly). Other mechanisms such as inclination-pumping during migration (e.g. iLee and Thommesl 
2009) may also contribute. 

In addition, previous dynamical studies have suggested that systems of prograde and retrograde planets are mor e 
stable than standard prograde-prograde cases (|Gavon and Boisl [20081 : iGavon et al.l 120091 : iSmith and Lissauerl 12009). 
essentially because the planets spend less time "close together" and thus perturbations are smaller. 

Given the observational evidence for highly inclined and retrograde systems, and the suggestion of enhanced stability 
in such systems, we investigate and quantify the hypothesis that retrograde systems will have a significantly reduced 
TTV profile compared to a standard prograde case. 

2. METHODOLOGY 

2.1. Integration and TTVs 

We consider a fiducial system consisting of two planets: a transiting Hot- Jupiter (with mass Mi = Mj np ^ er , 
eccentricity e\ = and semi-major axis a\ = 0.05 AU) and a small non-transiting, Earth-mass planet which perturbs 
the transit times of the larger planet. We evolve each system forward in time for ~ 10 years, equating to ~ 1,000 
transits of the Jupiter-mass planet (assuming that Hot- Jupiter remains transiting throughout). 

We perform the n-body integrations using a conservative Bulirsch-Stoer integrator, derived from that of Mercury 
(Chambers 1999). We use a barycentric coordinate system and limit the time steps to no more than 0.04 times the 
orbital period. After each time step, we test whether the projected star-planet separation (A) passed through a local 
minimum and the planet in question is closer to the observer than the star. Each time these conditions are met, we 
find the nearby time that minimizes A via Newton- Raphson iteration and increment an index i. If the minimum A is 
less than the stellar radius, then we record the mid-time of the transit, U. Observations of transit times are perturbed 
by the light travel time, S Uu(i) ~ — (r p (U) • r^ os )/c, where v p (ti) is the barycentric vector of the planet at time 
ri os is a unit vector pointing to the observer, and c is the speed of light. The observable transit time variations are 
calculated as S t(i) = t{ + S tm{i) — i x P — to, where the constant P and to are determined by linear least squares 
minimization of t(i)) 2 . We neglect any motion of the stellar center between the time of light emission and the 

time of transit. 

2.2. TTV Maps 

In £13.21 and 13.31 we produce TTV contour maps in (22,^2) and (^,^2) parameter-space. These plots are produced 
by fixing the mass and the initial orbit of the transiting hot- Jupiter (ei = & a\ — 0.05), and then varying the orbit 
of the Earth-mass perturber. 

For each set of configurations - i.e. semi-major axis, <22, eccentricity, e2, & inclination, %2 of the perturbing planet 
- we conduct 5 simulations randomising the other orbital elements (argument of pericentre, cj, longitude of ascending 
node, f2, & initial mean anomaly M). 

We then calculate the RMS of St(i) for each a, e, i configuration. The required contour plots are then plotted using 
an approximately logarithmic color scale (shown in Figs. [2] and [3]). 

We restrict o ur simulations to the region defined by the prograde coplanar three-body Hill stability criterion of 
iGladmanl (fl993h . This ensures that all of the i = simulations a re stable, and ensure s for t he i > systems that 
we have covered all of the guaranteed stability zone and more (see IVeras and Armitagd (j2004f ) for a discussion of the 
stability limit as a function of inclination). As an additional check on systems with extremely large TTVs close to the 
Hill-stability criterion boundary, we conducted 1,000 longer-term Mercury simulations. In most cases the Earth-mass 
planet experiences significant changes in semi- major axis within 1,000 years, due to the "kicks" at pericentre passage. 
Thus, although these planets do not collide, the less massive planet can experience some orbit variation. 

3. RESULTS 

3.1. Example TTVs For a Hot- Jupiter System 

In Fig QJi we show a circular hot- Jupiter at 0.05 au being perturbed by an Earth-mass companion with e2 = 0.02 
located on the external 2:1 mean- motion resonance (MMR). We find that although the prograde systems (0 < i < 90 °) 
have TTV profiles that differ significantly, their overall amplitudes are similar. In contrast, the amplitude of the 
i = 180° signal is very low (~ 1 second). The results for the 3:1 case in Fig[lf> are qualitatively similar. 

In contrast, for the 5:1 MMR system (Fig[lfc), while the amplitude of the i = 180° signal is extremely low, there 
are now significant differences in the behaviour at smaller angles: the i = and i = 180 now have similarly low 
amplitudes, with an approximate maximum occurring at i = 90°. 

Away from resonance (Figs [Dl, []J & [if), we again find that the i = 180° signal is generally significantly lower 
than that observed for cases with i < 90°, but that the detailed dependence on inclination is far from obvious, with 
some cases (Figs [T]i & [if) having the highest amplitude at i = 0°, while the low eccentricity plot of Fig. has a 
particularly low amplitude signal at i = . 

We note that there are a number of frequencies evident in FigHJ with periods from a few months to >> 10 years. 
In general, the short-term oscillations tend to be driven by variations in e2, leading to variations in the distance of 
closest approach between the two planets. The longer term quadratic trends are primarily due to outward drifts in 
a<ii 1 a\ as a result of kicks at close approach. 

The large number of different frequencies evident in the plots also serves to highlight the importance of considering 
the sampling period over which observations and simulations are to be conducted and compared. 
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a) 2:1 Resonance, e=0.025, i=0 
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Fig. 1. — The effect of varying inclination at fixed semi-major axes and eccentricities. Prograde orbits with: i = 0° have solid red lines; 
i = 45 ° have dotted green lines; 90 ° have dashed blue lines; while retrograde i = 180 ° orbits have dot-dashed black lines. In (a) we 
show the TTVs for a system on/near the 2:1 MMR (a 2 /ai = 1.59,P 2 /Pi = 2.0, e 2 = 0.02), In (b) the system is on/near the 3:1 MMR 
(a 2 /ai = 2.08,P 2 /Pi = 3.0, e 2 = 0.12), In (c) the system is on/near the 5:1 MMR (a 2 /ai = 2.92, P 2 /Pi = 5.0, e 2 = 0.05), In (d) we the 
system is just outside the 2:1 MMR (a 2 /ai = 1.63, P 2 /Pi = 2.08, e 2 = 0.05), In (e) we show a a non-resonant system at low eccentricity 
(a 2 /ai = 3.6, P2/P1 = 6.83, e 2 = 0.15), In (f) we show a non-resonant system at high eccentricity (a 2 /ai = 3.6,P 2 /Pi = 6.83, e 2 = 0.55). 
The results are complex, with different systems exhibiting a number of differing frequency contributions, with many contributions having 
different dependencies on inclination. However it is clear from these plots that in all the examples shown, the TTV signal for the i = 180 
case tends to be significantly lower in amplitude than for other inclinations. 

Furthermore, TTV profiles may change in amplitude and frequency by seve ral orders of magnitude simply due to 
a variation in the initial values of the mean anomalies ([Veras and Fordll2QQ9f ). This means that while approximate 
expressions for TTV amplitude dependencies on mass, semi-major axis separation, etc, can serve as a useful guide, 
individual system variations can be far removed from the mean. We do not intend this statement as a criticism of 
previous work, but rather as a cautionary note against over-reliance on approximations when analyzing individual 
systems which can exhibit great sensitivity to initial conditions. We do not seek in this paper to address in detail all 
such issues, but instead defer a more detailed consideration to a companion paper Veras et al. 2010, in prep., where a 
more thorough examination of the coplanar prograde case will be presented, allowing a much more detailed exposition 
of the huge number of variables which can affect the TTV signal. 



3.2. Eccentricity- Inclination Contour Plots 

To better understand the TTV amplitudes as a function of inclination, for each of the semi-major axis plots in Fig [I] 
we perform a search of the (e2, 12) parameter space, producing contour plots in Fig [2] showing the RMS TTV amplitude 
variation at that semi-major axis separation as the eccentricity and inclination of the outer planet is varied. 

In the 2:1 case of Fig [2^ the resulting contour plot is rather uniform, with the only area of significant reduction 
being in the e2 < 0.1, i ~ 180° region. An approximately similar result can be seen for the 3:1 case. 

For the 5:1 case (FigEfc), the situation changes: the region of low amplitude close to i = 180 significantly expands, 
while another low amplitude region emerges at e^ < 0.1 & i < 45 °. 

For Fig[2]i the perturber is located very close to (but not on) the 2:1 MMR. The plotted behavior is particularly 
rich, showing huge variations in the predicted TTV amplitude as functions of both eccentricity and inclination. 

Finally, the non-resonant case in Fig [2^ with 0,2/0,1 = 3.6, shows a fairly simple decline in TTV amplitude as a 
function of inclination. 

3.3. TTV Maps for Varying ^ 

We now look at a different projection of the data, looking in the (^-, e2) plane at the TTV structure as a function 
of inclination. 

In Fig [3^i we reproduce the coplanar prograde results of lAgol et al.l (|2005f ). showing the pronounced "flames of reso- 
nance" (jVeras and Fo rd 2009). The majority of the non-resonant parameter-space has TTV signals of comparatively 
low amplitude (< 30 seconds). It is only when eccentricities become very high, or the planets are in / near MMRs, 
that the TTV signal amplitudes become large. 

As the relative inclination of the planetary orbits is increased, for low inclinations changes in TTV amplitude are 
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Relative Inclination of the Outer Body (12)- 

Fig. 2. — Plots at fixed (02/01) showing the median RMS TTV amplitude as a function of e2 & 12. The contour key is displayed in the 
bottom right-hand corner. The semi-major axes selected for plots (a) - (e) are the same as those of Fig[TJi - Fig[T^. We can see significant 
differences in behaviour as a function of inclination depending on (i) position on / off resonance, (ii) the type of any resonance, and (iii) 
the eccentricity of the the perturber. 

relatively subtle; for the z = 30°&z = 45° cases (Fig [3t & [3]i) the RMS TTV signal is reduced in some high- 
eccentricity locations, corresponding to the reduction previously seen in the close-to-resonance map of Fig EJi. In 
addition, we can also see that spikes around higher resonances st art t o drop lower down, i.e. the low eccentricity 
amplitude increases as the inclination decreases, an effect noted in §3.21 for the 5:1 resonance plot of Fig[2fc. 

For larger relative inclinations the expected TTV signal amplitude starts to drop across much of the (^^2) param- 
eter space (see the i = 90 plot in Fig[3]i). The reduction in high amplitude signals at high eccentricities is the most 
marked decrease, but at any given eccentricity one tends to find a reduction in amplitude (shifting of the contours) 
compared to the i = plot. 

When the systems are pushed into retrograde orbits, the difference becomes more pronounced (E.g. i = 135 in Fig 
[3k and i = 180 in Fig[3]F). For these inclinations we find a reduction in amplitude across large swathes of the plot, 
with significant amplitude remaining only near MMRs. In addition, the TTV amplitude and the width of regions with 
large TTVs (near MMRs) is markedly reduced. 

3.3.1. External 2:1 MMR 

Since highly inclined systems only have large RMS TTVs in / near MMRs (Fig[3f), we examine in more detail 
one such resonant region. The external 2:1 MMR is known to be of importance for exo-planet systems from both 
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Semi-major Axis Ratio ((12/ai). 

Fig. 3. — Median RMS TTV amplitudes in the (02/01,62) plane for various 12- The contour key is shown at the top-left in plot (a). In 
each (02/01,62) plot, we see a strong increase in amplitude (i) at high eccentricity, and (ii) close to MMRs. As the inclination is increased 
there is an obvious general decrease in TTV amplitude, particularly pronounced for the retrograde orbits. 



theoretical (e.g. Sand or et al.ll2QQ7h and observational (e.g. iTinney et al.l 1200(1 ) studies. 

On and adjacent to the 2:1 MMR region, the TTV signal is extremely sensitive to changes in semi-major axis and 
/ or eccentricity. In the co-planar prograde case (Fig UK) j a fractional increase of just a few percent in either a 2 or 62 
causes the TTV amplitude to jump by more than an order of magnitude. 

As the inclination increases to 30 and 45 °, we find that the central region of the MMR maintains an approximately 
constant TTV amplitude. The regions just outside the resonance feel the strongest effects. E.g. regions near 0,2/0,1 = 
1.54 show TTV amplitudes reduced by an order of magnitude for many eccentricities (compare Fig UK and Fig Eh). 
This decreases the width of the region of parameter space around the 2:1 MMR resonance with TTV signals > 30 
seconds. In addition, the general reduction in "background amplitude" causes the minor resonances to become more 
prominent: E.g. 02/01 ~ 1.69. 

As the inclination increases to 90° and 135°, TTV signals > 30 seconds are restricted purely to the resonance 
regions, with regions away from resonance having signals < 30 seconds almost everywhere, even at high eccentricity. 
As the inclination reaches 180°, and the planets are now again coplanar but with completely counter-rotating orbits, 
we find TTV amplitude over 100 seconds only in very narrow regions around the first- and second-order MMRs. 

Finally, we note that as we move from Fig 0k to Fig|4f, the amplitude in regions slightly offset from resonance appears 
to increase slightly at high eccentricities. However this only occurs at low order resonances (small semi-major axis 
separations - see Fig[3k and Fig[3F) as well as Fig El More generally, away from such regions (see Fig[3]F), retrograde 
orbits tend to produce the smallest TTVs. 
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Fig. 4. — Close-up of the 2:1 MMR region, demonstrating extreme sensitivity to variations in 02/01, e2 Sz %2- The contour key is the 
same as that used in Figs. 

3.3.2. Internal Perturbers 

Finally, for completeness we consider an Earth-mass planet on an orbit interior to the Jupiter-mass planet, keeping 
the more massive planet at a\ = 0.05AU and e\ = 0, and then monitor the TTVs of the massive planet due to the 
interior perturber. 

From Fig [5] it seems that the background (non-resonant) regions now retain an approximately constant TTV signal 
irrespective of inclination, while there is a tendency for the TTV signal to become stronger around the MMR regions 
as the inclination is increased towards 90 , but then drop away to almost nothing as the inclination further increases 
towards i = 180°. 

4. CONCLUSIONS 

We have investigated the TTV signals for systems of highly inclined and retrograde planets. We find that: 

• In the vicinity of exterior MMRs the inclination dependence is complex: low order resonances maintain a high 
TTV amplitude for all regions i < 170 °, declining in amplitude only for low eccentricity cases close to i = 180 °, 
where-as higher order resonances display an increase in TTV amplitude as inclinations rise from to 45 degrees. 
Moreover, the regions immediately adjacent to MMRs show extreme sensitivity to changes in perturber a, e and 
i. 



Exterior perturbers away from resonances tend to show a slow decrease in TTV amplitude with increasing 
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Fig. 5. — Internal perturber parameter space. Top; a general view. Bottom; a close-up of the 1:2 MMR. The contour key is the same 
as in Fig@] We find that the behaviour in the internal region differs from the external region: the background stays constant and there is 
now a strengthening of the resonant regions as the inclination increases towards 90 . 

inclination, although regions adjacent to resonances can show remarkably complex behaviour. 

• Perturbing planets on interior orbits display a slightly different behaviour: Away from resonance the amplitude 
remains approximately constant with inclination, but around MMRs the perturbations become stronger as the 
inclination increases towards 90° before decreasing again beyond 90°. 

We note that: 

• Absence of evidence is not evidence of absence: Planets in retrograde orbits should be expected to produce 
markedly reduced TTV signals as compared to the standard prograde case. For an Earth-mass perturber in an 
anti- aligned orbit, almost the entirety of the sample parameter space would result in a very small or undetectable 
TTV signal (unless the planet happened to be fortuitously located precisely on a strong MMR). 

• Retrograde orbits may be a natural way to explain transiting systems in which little or no TTV signal is observed, 
but in which the radial velocity observations point towards the existence of an additional planetary companions 
(E.g. GJ 436, HAT-P-13, etc). 

• In addition to the TTV considerations in this work, it is important to acknowledge that inclined orbits in 
multi-planet systems will precess, leading to variations in transit duration. More work is required to try and 
understand whether a combination of TTV signals with transit duration variation signals could remove some of 
the degeneracies inherent in this problem. 

• Many additional dependencies (E.g. sampling period, perturber mass, relative mean anomaly, etc, etc) can signif- 
icantly alter the expected TTV signal for a particular system. We defer the provision of a detailed investigation 
of such matters to a companion paper, Veras et al. 2010, in prep. 

• Future work could also investigate in more detail the regions above the Hill stability curve, seeking to identify 
the (likely large) TTVs for any stable systems which exist in that region. 
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